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Аннотация 

Введение. На сегодняшний день существует множество способов амплификации нуклеиновых кислот и у 
каждого способа есть ряд достоинств и недостатков. Одним из наиболее популярных способов является петлевая 
изотермическая амплификация (Г.оор-те1ае4 150 'егта! АМРИйсайоп, ГАМР). В отличие от термоциклических 
реакций, таких как ПЦР (полимеразная цепная реакция), для которых требуется смена трех температурных 
режимов и дорогостоящее оборудование, в Г.АМР вся реакция проходит при одной температуре и с максимальной 
на данный момент скоростью. Важным компонентом проведения ГАМР-амплификации являются праймеры 
(обычно 20-25 нуклеотидов), которые необходимо подбирать к определенному участку нуклеотидной 
последовательности. Известно, что последовательность ДНК содержит четыре нуклеотида: А — аденин и 
Т — тимин, Г — гуанин и Ц — цитозин. Вариантов перестановок этих нуклеотидов огромное множество, и 
проанализировать вручную такое большое количество данных практически невозможно, поэтому возникает 
необходимость в использовании современных компьютерных технологий. Для дизайна ПЦР-праймеров 
предложено более 150 компьютерных программ, в то время как для ГАМР-праймеров их менее 10, и каждая из 
них имеет ряд недостатков, например, по длине анализируемого участка. Поэтому целью данной работы является 
разработка новой отечественной компьютерной программы дизайна специфичных праймеров именно для ГАМР. 
Материалы и методы. В основе алгоритма поиска праймеров лежит линейный поиск подстроки в строке с 
учетом критериев подбора праймеров для ГАМР. Программный комплекс дизайна ГАМР-праймеров разработан 
на языке программирования Ру#фоп. Для работы с различными ДНК и РНК использовалась библиотека оРуоп, 
а для разработки интерфейса — фреймворк ©. 

Результаты исследования. Предложена модификация метода прямого перебора с использованием трафаретного 
подхода, учитывающего @С-состав и температуру отжига праймеров в зависимости от их структуры. Разработан 
комплекс программ с дружелюбным интерфейсом, учитывающий критерии дизайна праймеров: получены 
свидетельства о регистрации программ для ЭВМ (ГАМРитег$ 10 № 2022617417 от 20 апреля 2022 года, 
ГАМРитегз 10 1оор № 2023662840 от 14 июня 2023 года). Программа есть в открытом доступе по адресу: 
В рь://о баб. сот/Везу/.АМРгитегз 1 

Обсуждение и заключение. Разработанные программные комплексы могут использоваться для исследований и 
анализа в области молекулярной биологии и генетики, для создания диагностических тест-систем, 
обеспечивающих высокую чувствительность и достоверность обнаружения специфических ДНК и РНК. 
Программные комплексы могут применяться в научно-исследовательских институтах и лабораториях, 
занимающихся амплификацией нуклеиновых кислот. Результаты оценки подобранных наборов праймеров для 
реакции ГАМР апробированы, и эффективность рабочих наборов с помощью программы ГАМРгитегз 19 
доказана экспериментально на примере обнаружения генетического материала коронавируса ЗАВ $-Со\У-2. 
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Введение. Амплификация нуклеиновых кислот является ценным молекулярным инструментом не только в 
фундаментальных исследованиях, но и в прикладных областях, например таких, как диагностика инфекционных 
заболеваний, наследственных патологий, установление родства и т.д. В настоящее время методы амплификации 
интенсивно развиваются, расширяются и области их применения. Самым популярным и наиболее часто 
используемым методом амплификации является полимеразная цепная реакция (ПЦР) [1]. ПЦР — это реакция, 
проходящая при трех различных температурных режимах: денатурация (95 °С), отжиг праймеров (от 50° до 60 °С), 
элонгация (72 °С). Для быстрой смены этих режимов необходим специальный прибор — ДНК-термоциклер [2]. При 
этом смены температурных режимов в амплификаторе происходят не моментально, а начинаются лишь по 
достижении нужной температуры, и это ведет к искусственному сдерживанию протекания реакции. Как правило, 
продолжительность ПЦР составляет 1-1,5 часа. 

Вторым по популярности методом амплификации является петлевая изотермическая амплификация (Гоор 
АМРИйсаноп, ГАМР) [3]. Для проведения ГАМР достаточно водяной бани или термостата, так как реакция 
проходит при одной температуре, и первые результаты можно увидеть уже через 15 минут. 

Как для ГАМР, так и для любого другого вида амплификации ключевым компонентом являются праймеры, 
которые представляют собой короткую последовательность нуклеиновой кислоты и служат отправной точкой 
для увеличения копий определенного участка ДНК. Именно праймеры определяют, какая последовательность 
ДНК будет копироваться. 

Основным отличием проведения ГАМР является количество праймеров. Для обычной ГАМР необходимо 
минимум четыре праймера (два внешних, два внутренних), в то время как для обычной ПЦР достаточно двух 
(прямой, обратный). 

Для повышения специфичности и точности реакции важно правильно подобрать праймеры. Для 
автоматического подбора праймеров для ПЦР разработано более 150 различных компьютерных программ, 
которые позволяют подбирать праймеры для любых модификаций этой реакции [4]. Однако для ГАМР таких 
программ крайне мало, не более десяти, и всего две из них доступны онлайн. Также эти программы имеют ряд 
недостатков, таких как ограничения по длине анализируемой последовательности, не исключают возможности 
образования гомо- и гетеродимеров праймеров, повторов нуклеотидов в одном праймере. И ни одна программа 
не учитывает сближенное расположение праймеров в одном наборе, что, в свою очередь, снижает качество 
праймеров и точность результатов реакции [5]. 

Вследствие этого актуальной задачей является разработка новой компьютерной программы, позволяющей 
подбирать (моделировать) качественные наборы праймеров для ГАМР-амплификации с ужесточенными 
условиями подбора праймеров для нуклеотидных последовательностей любой длины. 

Материалы и методы. Авторы оригинальной статьи [3] предложили использовать два внешних, 
ЕЗ (Еог\ага), ВЗ (ВасК\ага), и два внутренних праймера, ЕТР (Еог\уаг4 Гапег Ргитег), ВТР (ВасК\г’агА Гппег Ргипег). 
Предполагалось, что внутренние праймеры имеют удвоенную длину (ЕР: Е1с/Е2, ВПР: В1с/В2) и отжигаются на 
четырех участках нуклеотидной последовательности. Схематично расположение ГАМР-праймеров можно 
увидеть на рис. 1. 

Внешние праймеры необходимы только на начальном этапе. Они предназначены для того, чтобы ограничить 
анализируемый участок нуклеотидной последовательности и образовать одноцепочечную структуру этого 
участка. Пара внутренних праймеров, Е1с и В1с, начинают свою работу уже на втором этапе, поскольку они 
отжигаются после образования новых цепей ДНК. 
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Рис. 1 Схематичное расположение мест отжига внешних и внутренних праймеров 
для ГАМР-амплификации [6] 


Позже теми же авторами был предложен модифицированный метод, предлагающий использование не 
четырех, а шести праймеров, отжигающихся уже на восьми участках целевой нуклеотидной 
последовательности [7]. Предлагалось добавить еще два петлевых праймера (Гоор В, Гоор Е), которые должны 
вступить в реакцию на третьем этапе после образования гантелеобразной структуры ДНК и отжигаться между 
участками Е1/Е2 и В1/В2 соответственно. Использование дополнительных праймеров подразумевает увеличение 
чувствительности и достоверности реакции. 

Любая реакция амплификации имеет свой порог чувствительности, и разброс данного показателя весьма 
велик ввиду того, что на протекание как ПЦР, так и ГАМР влияет множество факторов. В некоторых работах 
отмечается, что ГАМР по чувствительности значительно превосходит ПЦР. Например, авторы работы [8] 
утверждают, что ГАМР чувствительнее ПЦР в 10 раз. Авторы других работ выявили, что ГАМР чувствительнее 
ПЦР 100 раз [9], ау некоторых этот показатель достигает 1 000 раз [10]. 

Помимо чувствительности, у любой реакции амплификации есть еще один не менее важный показатель — ее 
специфичность. И здесь к ГАМР-амплификации в последнее время стали появляться вопросы [11], в том числе 
из-за возникновения так называемых праймерных гомо- и гетеродимеров, исключить которые в этой реакции 
сложнее, чем в ПЦР, ввиду большего числа используемых праймеров и их увеличенной длины [12]. 

Для успешного проведения амплификации необходимо правильно подобрать праймеры. При использовании 
метода ГАМР основная сложность заключается в моделировании праймеров с учетом всех рекомендуемых 
условий, а именно: 

1) длина праймера (18-35 нуклеотидов для внешних праймеров, 30-55 нуклеотидов для внутренних); 

2) содержание гуанина (С) и цитозина (С) (ОС-состав в пределах от 40 до 60 %); 

3) оптимальная температура отжига праймеров (55—65 °С); 

4) сближенное расположение праймеров в одном наборе: средний размер ампликона (120—220 п. н.); 

5) исключение образования димеров праймеров; 

6) исключение повторов нуклеотидов в одном праймере (не более трех). 

В таблице 1 приведены краткие характеристики и возможности наиболее популярных программ дизайна 
ГАМР-праймеров [5]. 
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Таблица 1 
Краткие характеристики популярных компьютерных программ дизайна 


ГАМР-праймеров 
= |) 
ы [ай я 
5 г д р Е Ра Яо ь- 
5 . Е : Е 
Е РНЕ |218 ВНЕ | Е 
я зано боя ло = 5 = ан 
ь: Е: ВИ> С-В БЕ о © © 
5. вне ва ве он 
Е Е > вв Ее ЗЕ Ф 
ы в. -® З - 
м | 
Ргипег Ехр|огег (\У4, \У5) 
японской компании Е\Жеп СБеписа1 Тауа до 2 000 Е нет да да 
Со. ГТО, Токио 
ЕаяРСК 
. . 12 — 
финской фирмы Ргитег О1еИа| Та о 988 Нот к __ 
ОТАРО С, 
Шанхайского СОРА С, Ре! не ограничено во ео нет да 
университета Цзяо Тун, Китай 0$ Глпих 
ГАМР Пез1 
и о . | до 15 000 000 + только Тм да нет 
американской фирмы Ргепиег В1озой 
Гатргип Руфоп — — — __ 
доступа 
МЕВ ГАМР, 
№ № ем Епапа В1ю]аб$, Ипсвиче, штат Тауа от 100 до 2 000 + — да да 
Массачусетс, США 
ГАМРитег 10. 
ГАМРитег 10 -Т0ор, Руоп 3.10 | не ограничено + да да да 
Уфа, Россия 


Дизайн праймеров для ГАМР-амплификации представляет собой весьма трудную задачу и требует 
разработки специальной компьютерной программы с должным функционалом, учетом всех рекомендуемых 
условий и с возможностью расширенного подбора праймеров и удобным интерфейсом. 

Компьютерная программа дизайна праймеров разработана на языке программирования Ру®оп. Данный язык 
имеет библиотеку МоРу®оп, которая позволяет работать с нуклеотидными последовательностями, а также 
фреймворк О! для разработки интерфейса. 

Результаты исследования. С учетом особенности строения нуклеотидных последовательностей и критериев 
подбора ГАМР-праймеров была предложена модификация метода прямого перебора с использованием 
трафаретного подхода, учитывающего СС-состав, температуру отжига праймеров и позволяющего снизить 
сложность перебора. 

Как известно, @С-состав праймеров должен быть в пределах от 40 до 60 %. Это один из важных критериев 
подбора ГАМР-праймеров, который зависит от самой анализируемой последовательности, длины праймера и 
частично влияет на температуру отжига (Тт, °С). 

В данной работе для расчета температуры отжига праймеров используется формула: 


Тт = 81,5+16,6. (ов [ Ма" |)+0,41-(%6+%С)-— 5481, (1) 
где [№ — молярная концентрация ионов натрия; (%С+%С) — СС-состав в анализируемой 
последовательности, выражается в процентах; С — длина праймера. За ее основу была взята известная 
зависимость: 

Тт = 81,5+16,6- [108 [№ ])+0,41.(%6+%С)- 600Леней. (2)! 


Формула (1) была подобрана эмпирически, расчетные значения сопоставлялись со значениями, получаемыми 
с помощью удобной утилиты ОНгоАпа|утег?, которая предоставляет качественный подбор праймеров для любых 


1 Ойго Сис: ОПзописеонае Ргорегиеу Сайсшаиог. ОВТ:: Вр://ЛБлоюо15 пис. погуееги.еди/ОНзоСас.Вии! (дата обращения: 10.12.2023). 
2 ОПооАпаугегтм Тоо]. ОВТ: Вирз:/еи1А па. сот/расез/оо15/оНеоапайухег гепипи =%2Еса!с%2Еапайутег (дата обращения: 10.12.2023). 
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видов амплификации нуклеотидных кислот. Определив длину праймеров, СС-состав и температуру отжига, 
можно найти все возможные праймеры в нуклеотидной последовательности любой длины, которые будут 
удовлетворять заданным критериям. 

Если представить праймер как подстроку, а анализируемую нуклеотидную последовательность как более 
длинную строку, данную задачу можно представить как перебор всех возможных вариантов (прямой поиск), но 
усложнив его за счет расчета СС-состава и температуры отжига праймеров. 

На рис. 2 представлена полная блок-схема работы алгоритма прямого перебора, учитывающая: 

— длину праймеров (п. о.); 

— ОС-состав, %; 

— температуру отжига праймеров, Тт, °С; 

— гомодимеры на обеих цепочках ДНК. 

Общая сложность модифицированного алгоритма в худшем случае: О(т-п), где п — длина праймера; 
т — длина нуклеотидной последовательности. Необходимо понимать, что время работы алгоритма напрямую 
зависит от того, насколько часто встречаются удовлетворяющие требованиям нуклеотидные фрагменты. 

В таблице 2 приведены данные поиска всех возможных праймеров в нуклеотидных последовательностях 
различной структуры по времени. 

Далее из всех найденных праймеров необходимо сформировать наборы, учитывая сближенное расстояние между 
праймерами, гетеродимерность в одном наборе, а также минимальную разницу температур отжига праймеров. 

Схема дизайна формирования праймеров в наборы для ГАМР представлена на рис. 3. 


Цикл поиска праймеров 
Анализ первых и - 1 нуклеотидов 
СС, расчёт флагов 


В пределах ь 
Возврат массива праймеров 
диапазона поиска 


Расчёт флагов и @С для 


нуклеотида текущего индекса 


М Расчёт параметров 


нет да 
Е Проверка на гомодимеры [. 


нет да 
Нет гомодимеров? 
Добавление праймера 


Проверка на гомодимеры 
с комплементарной цепи 


нет да 
Нет гомодимеров? 
Добавление праймера 


Расчёт ОС и флагов для 


первого и последнего 
нуклеотида 


Рис. 2. Полная блок-схема работы алгоритма 
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Таблица 2 
Время, затраченное на поиск подходящих праймеров в нуклеотидных 
последовательностях разной длины 
Название генома НО Поиск праймеров, с 
последовательности, п. о. 
ЗАК5-Со\-2 29 903 0,31 
Езспенсва ушиз Т4 168 903 1,73 
Мусор!азта 580 076 5,43 
Нейсобацег руют 1624 458 18,11 
ЕзсБепсва сой 4 641 652 71,68 (1,2 мин) 
Саепограб 411$ 100 286 401 1082,53 (18 мин) 
Начало 


нет да 


Загрузка из файла? 
Загрузка файла 
Ввод 
последовательности 
Выбор хромосомы 


Установка конфигурации 
Поиск праймеров, 
расчет длины, Тт, О@С-соства 
Сортировка праймеров 
в наборы 
Вывод наборов 
на экран 

Конец 


Рис. 3. Блок-схема формирования наборов праймеров для ГАМР 


На входе программа считывает нуклеотидную последовательность, либо загружается нужный файл, либо 
вставляется фрагмент через буфер обмена. Далее устанавливаются конфигурации, такие как длина праймеров, 
СС-состав, температура отжига праймеров и разница температур в одном наборе. Затем происходит поиск всех 
возможных праймеров, удовлетворяющих установленным конфигурациям, подобранные праймеры сортируются 
в наборы и выводятся на экран пользователю. 

Работа компьютерной программы: 

1. Загрузка файла (простой текстовый формат, формат ЕАЗТА, СепВапК) или фрагмент последовательности 
через буфер обмена. 

2. Поиск всех возможных праймеров: праймеры объединяются по следующим критериям: 

— длина анализируемого участка; 

— расстояния между праймерами (ЕЗ/Е2 — 1-10 нуклеотидов, Е2/Е1с — 10-25 нуклеотидов, 
Е1с/В1с — 0-30 нуклеотидов); 

— разница температур отжига праймеров (33); 

— гетеродимерность (Беегодитег). 

В случае соблюдения всех вышеперечисленных условий набор считается рабочим. 

3. Вывод смоделированных наборов праймеров на экран пользователю и/или сохранение в файл. 
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Программа для ГАМР-амплификации зарегистрирована в Реестре программ для ЭВМ под названием 
ГАМРитегз 10, № 2022617417 20 апреля 2022 года и ГАМРитегз 10-юор, № 2023662840 14 июня 2023 года. 
Код программы в открытом доступе. 

Разработанный программный продукт обладает дружественным и интуитивно понятным интерфейсом, что 
позволяет использовать его непосредственно конечным пользователям — экспериментаторам, занимающимся 
ГАМР-амплификацией. 

Обсуждение и заключение. Количество выданных наборов праймеров зависит от заданных параметров 
поиска. Чем строже параметры, тем меньше наборов будет найдено. В случае жестких ограничений программа 
может не выдать ни одного набора. Количество наборов праймеров при различных параметрах подбора для 
генома бактериофага лямбда [13], длина которого составляет -48 500 нуклеотидов, представлено в таблице 3. 

Для относительно коротких нуклеотидных последовательностей (до 2 000 нуклеотидов) подбор праймеров 
занимает менее секунды. С увеличением длины последовательности длительность поиска праймеров 
увеличивается экспоненциально. 

Таблица 3 
Количество наборов праймеров для генома бактериофага лямбда в зависимости 
от заданных параметров отбора 


сс, % АТ, °С ине длина мтифишируомкио о Н. 

5 213 119 3 
а 2 198 184 3 

5 132 80 0 
+82 2 116 64 0 

5 195 185 4 
о р) 181 164 4 

5 160 147 4 
а р) 134 105 0 


На рис. 4 представлено влияние длины нуклеотидной последовательности на длительность подбора наборов 
праймеров. Приведены данные для нуклеотидной последовательности бактериофага лямбда на ноутбуке с 
параметрами: процессор Пие(В) Соге(ТМ) 17-10750Н СРО, 2.60СН?, 6 ядер. 16 ГБ ОЗУ. Указаны параметры 
мягкого подбора праймеров (40—60 % СС, АТ = 5, длина анализируемого участка — до 300 п. н.). 


3000 


2000 


1000 


Продолжительность поиска наборов, с 


1к 2к 5к 10к 20к 40к 


Длина ну клеотидной последовательности,. н 


Рис. 4. Влияние длины нуклеотидной последовательности на длительность подбора наборов праймеров 


3 ГАМРитег$-Ю. ОВТ.: Вирз://оаЬ.сот/Везву/Г.АМРгипетз-{О/Боб/таЛатр/ ат 1атр.ру (дата обращения: 10.12.2023). 
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Следует отметить, что продолжительность поиска праймеров зависит от мощности компьютера. 

Для сравнения и определения качества смоделированных наборов праймеров был проведен ряд натурных 
экспериментов по обнаружению РНК коронавируса ЗАВ$-СоУ-2, длина которого - 30 000 нуклеотидов. Для 
этого к одному и тому же участку нуклеотидной последовательности коронавируса были подобраны наборы 
ГАМР-праймеров с помощью программы ГАМРитег$ 10 и двух популярных и доступных онлайн-утилитов от 
компании Ме\ Еп1апа В1о1аб5 (МЕВ ГАМР)“ и РитегЕхр|огег°. Обозначения 1, М и Р соответствуют наборам 
праймеров ГАМРитегз 10, МЕВ ГАМР Ритег Оезеп и РитегЕхр|огег; «+» — образцы содержали РНК 
коронавируса ЗАВ 5-СоуУ-2, «—» — образцы не содержали нуклеиновых кислот. На рис. 5 приведены кривые 
данного эксперимента. 


О О О | 
А] 

1000 - 
ГАА 


Флуоресценция, отн.ед 


Продолжительно сть реакции, мин 


Рис. 5. График сравнительного эксперимента: 
Г, — праймеры, полученные с помощью ГАМРгипегз 14, № — МЕВ ГАМР Ритег Оез12п; 
Р — РишегЕхр!огег, «+» — образцы содержали РНК коронавируса ЗАВ $-СоуУ-2; 
«—» — контрольные образцы без матрицы (не содержали нуклеиновых кислот) 


Праймеры, полученные с помощью РитегЕхр|огег, показали самый поздний подъем кривых 
амплификации (Р-+), по сравнению с МЕВ ГАМР (М№+) и ГАМРитегз 10 (Г.+). Праймеры, полученные с помощью 
ГАМРитегз 10, обеспечивали более поздний подъем кривых амплификации ([.+) по сравнению с (№+). Однако 
образцы, не содержащие РНК вируса (Р-—), выдали более поздний подъем по сравнению с набором (№-), в то время 
как (Г—) не показали никаких подъёмов даже после 50 минут, тем самым обеспечивая самую высокую надежность 
обнаружения вирусной РНК. 

Проведенные эксперименты показали более высокую точность и специфичность наборов праймеров, 
подобранных с помощью компьютерной программы ГАМРиштегз 10, вызванную снижением скорости 
протекания реакции с образцами отрицательного контроля. Кривые амплификации имели более поздние 
подъемы или не имели их вообще, даже по истечении 50 минут продолжительности реакции. 
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